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I. INTRODUCTION 

The calculation of the band structure of metals and 
insulators is a central task in solid-state theory. A com- 
monly accepted method for this purpose is the density- 
functional theory (DFT) which provides surprisingly ac- 
curate results for the band structure of many materi- 
als 1 . Furthermore, the DFT is an 'ab-initio' theory, i.e., 
it starts from the full Hamiltonian of a real system and 
does not require the introduction of any simplified mod- 
els. The only competing 'ab-initio' theory was Hartree- 
Fock theory, which has shown very many shortcomings as 
compared to DFT, such as gross overestimation of band 
widths and band gaps. 

However, from a theoretical point of view, the suc- 
cess of the DFT for band structures is rather astonish- 
ing because this theory is a generic approach to ground- 
state properties only. All results on energy bands are 
extracted from auxiliary one-particle dispersions which 
have no physical meaning at the outset. Indeed, some 
shortcomings of the DFT energy bands have become ev- 
ident very early, in particular the underestimation of 
the fundamental gap in semiconductors. In semiconduc- 
tors and insulators, the so-called GW approximation to 
the one-particle Green function has been put forward 2 . 
There, the single-particle self energy is calculated using 
a Green function based on the DFT wave functions and 
the screened Coulomb interaction. It turns out that the 
GW quasi-particle bands are more or less rigidly shifted 
against the DFT bands so that the band gap results of 
GW calculations for semiconductors and insulators agree 
much better with experiment. 

For materials with strong electron-electron interac- 
tions, the DFT results have not been too convincing, 
in particular for magnetic insulators and other strongly 



correlated electron systems. For the iron group met- 
als the discrepancies of DFT results to experimental 
data, e.g., angle-resolved photo-emission spectroscopy 
(ARPES), increase towards the end of the series, i.e., 
towards nickel; for a detailed discussion on the discrep- 
ancies between DF T results and experimental data on 
nickel, see Refs. l8l ll For the iron group metals, GW cal- 
culations did not yield significant improvements over the 
DFT results; for nickel, see Ref. U 

A proper description of solids with strong Coulomb 
interactions requires true many-particle approaches. In 
the past, the notorious difficulties of many-particle sys- 
tems have restricted such theories to the study of rather 
simplified model systems, e.g., the one-band Hubbard 
model. Therefore, a comparison with experiments on 
real materials could hardly be performed. Only recently, 
new non-perturbative many-particle methods have be- 
come available which have made possible the investiga- 
tion of more rea listi c many-particle models; see, for ex- 
ample, Refs. CTl ltlTl lsl In Ref. we introduced a class 
of Gutzwiller variational wave functions which allow us 
to study general multi-band Hubbard models. Expecta- 
tion values with these correlated electron states are eval- 
uated exactly in the limit of large spatial dimensions, 
D — ► co. When applied to nickel, the remaining min- 
imization problem is numerically non-trivial because of 
the large number of vari ational parameters; first results 
are reported in Refs. I'll 

The Gutzwiller variational theory provides an approxi- 
mate picture of the ground state but, in principle, it lacks 
any information about excited states. This drawback can 
be overcome in two ways. First, if we take for granted 
that the variational ground state is at least qualitatively 
close to the true ground state, we may use the varia- 
tional state as a starting point for the variational calcu- 
lation of excited states. In Ref. we have used this idea 
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to determine the spin-wave dispersion in ferromagnetic 
multi-band Hubbard models. We have successfully repro- 
duced the experimental observation that the low-energy 
spin excitations in itinerant ferromagnets are very simi- 
lar to those of a system with localized spins. Second, the 
calculation of the variational ground-state energy in the 
limit of infinite dimensions^ naturally leads to the defi- 
nition of an effective single-particle Hamiltonian which, 
in some limits, can also be derived in Slave-Boson mean- 
field theorjiiS. Very much in the spirit of the DFT we 
used the band structure of this effective Hamiltonian for 
a successful comparison with ARPES data for nickel. We 
have been able to resolve basically all of the LDA short- 
comings. 

Despite its success, the second approach still lacks a 
sound theoretical basis. In this work we derive the (vari- 
ational) quasi-particle dispersion referring back to Lan- 
dau's original ideas on Fermi liquids. The Gutzwillcr 
variational state is an illustrative example for a Fermi- 
liquid ground state: the Gutzwiller many-body correlator 
acts on the Fermi-gas ground state whereby energetically 
unfavorable configurations are gradually reduced. In the 
spirit of Fermi-liquid theory a quasi-particle excitation is 
readily viewed as a Gutzwiller-correlated single-particle 
excitation of the Fermi-gas ground state. The energy of 
this excitation is identical to the quasi-particle disper- 
sion in our original workS. Therefore, no revision of our 
previous numerical results on nickei&i is necessary. 

The evaluation of the variational energy in our method 
is exact only in the limit of infinite spatial dimensions. 
Our application to realistic three-dimensional systems re- 
quires that 1/D corrections are well controlled. As is 
known from the one-band model, these corrections are 
small for ground-state properties such as the (variational) 
energy or the effective mass of the quasi-particles at the 
Fermi surface. In this work we will present additional re- 
sults on 1/D corrections of the quasi-particle dispersion 
for the one-band Hubbard model. 

Our work is organized as follows. In Sec. [H] we sum- 
marize the basic ideas of Landau-Gutzwiller theory. In 
Sect. lIIll we discuss the variational ground state for multi- 
band Hubbard models. In Sec. IIVI we define Landau- 
Gutzwiller quasi-particles and derive their energy dis- 
persion. In Sect. [V] we calculate 1/D corrections for the 
quasi-particle dispersion of the one-band model. Our pre- 
sentation closes with short conclusions. Some technical 
details are deferred to the appendix. 

II. LANDAU-GUTZWILLER THEORY 

In second quantization the Hamilton operator for non- 
interacting electrons reads 

Hi = Cf ^<£j\<T' • y2<*?<;lr<-,:n ■ (1) 

Here, c£~ a creates an electron with combined spin-orbit 
index a — 1,...,2N (N = 5 for 3d electrons) at the 



lattice site i of a solid with L lattice sites. The electron 
density n = N c /L, is finite in the thermodynamic limit 
N e — > oo, L — > oo. 

For a translationally invariant system, as considered 
throughout this work, this single-particle Hamiltonian is 
readily diagonalized in momentum space. Its eigenstates 
are one-particle product states \*5f). In particular, the 
ground state \^q) is the filled Fermi sea where all single- 
particle states below the Fermi energy are occupied. All 
other eigenstates can be understood as particle-hole ex- 
citations of l^o)- 

One essential idea behind Landau's Fermi liquid the- 
ory is the assumption that the Fermi-gas picture remains 
valid qualitatively when electron-electron interactions are 
switched on; for an introduction, see, e.g., Ref. I 111 The 
Fermi-gas eigenstates transform adiabatically into those 
of the Fermi liquid while keeping their physical prop- 
erties. For example, the momentum distribution dis- 
plays a discontinuity at the Fermi energy both in the 
Fermi gas and in the Fermi liquid. Naturally, the prop- 
erties of the ground state and of the particle-hole exci- 
tations change quantitatively. Therefore, the excitations 
are called quasi-particles and quasi-holes in the Fermi 
liquid. 

Gutzwiller's variational theory closely follows the idea 
of an adiabatic continuity from the Fermi gas to the Fermi 
liquid. Let us introduce a general class of Gutzwiller- 
correlated wave functions via 

|# G ) = P G |$). (2) 

Here, |<E>) is any normalized one-particle product state. 
The Gutzwiller correlator 

p G = Uh G (3) 

i 

is a many-body operator which suppresses those configu- 
rations which are energetically unfavorable with respect 
to the electron-electron interaction. Therefore, the (ap- 
proximate) Fermi-liquid ground state 

|H/ G ) = P G |$ ) (4) 

evolves smoothly from a Fermi-gas ground state |<E>o) 
when the electron-electron interaction is switched on. In 
fact, this concept has been used by Vollhardlii^ to develop 
a microscopic theory for the ground-state properties of 
liquid 3 He on the basis of Gutzwiller's approach. 

In this work, we use Landau's idea to extend Gutzwil- 
ler's variational approach to quasi-particle excitations. 
In principle, this does not pose a big problem. Instead of 
a Fermi-liquid ground state \$o), we use single-particle 
excitations |<I>) of the Fermi gas in J2J) to define quasi- 
particle states. In Sect. IIVI we will give a proper math- 
ematical definition of a quasi-particle excitation. Here 
we point out that all restrictions of Fermi-liquid theory 
apply. For example, only the low-energy properties of 
metals, close to the Fermi energy, ought to be described 
in this way. Nevertheless, experiments on metals show 
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that well-defined but life-time broadened quasi-particle 
excitations can be found even for energies of about 10 eV 
below the Fermi energy. Therefore, the concept of quasi- 
particles and quasi-holes remains meaningful for those 
parts of the valence and conduction bands which are rel- 
evant in solid-state physics. 



III. VARIATIONAL ENERGY 

A. Multi-band Hubbard Hamiltonian and 
Gutzwiller variational states 

In the following we study multi-band Hubbard models 
where the electron-electron interaction is purely local, 



H = Hi 



E 



Hi,, 



(5) 



Here, the atomic Hamiltonian Hi- ai contains all possible 
Coulomb-interaction terms between electrons on site i, 



Hi ;a t — 



E 



2;(7l 1\<J2 «;CT3 Z;(T4 



(6) 



We assume that the eigenstates |r),; of the atomic Hamil- 
tonian have been determined 



Hi,, 



E^ 



r™ 1: r 



rrii-v 



(7) 



This is possible in all cases of interest, at least numeri- 
cally. In the following, the site index will often be sup- 
pressed as we are primarily interested in translationally 
invariant systems. 

The Gutzwiller theory allows us to study the Hamil- 
tonian J5j with an arbitrary number of orbitals^. In this 
work, however, we will restrict ourselves to the special 
case where non-degenerate orbitals belong to different 
representations of the respective point-symmetry group. 
For example, in cubic symmetry we allow for only one 
set of s, p, d(e g ) and dfag) orbitals. 

The Gutzwiller correlator 



is parameterized by 2 2N real numbers Ai : r- For an ener- 
getically unfavorable atomic configuration (Tj) the mini- 
mization will result in Ai ; r < 1 whereby its weight in |<J>) 
is reduced. 



B. Extrema of the variational energy 

In the limit of large spatial dimensions, the expectation 
value of the Hamiltonian JSJ in the wave function can 
be expressed in terms of the one-particle product wave 
function |$) and the expectation values 



(m r )* c 



(^ G |m r |^ G ) 
(*g|*g> 



(9) 



for all details, see Ref.0 After a lengthy calculation one 
obtains the following variational energy 



S a ,<r>(k) 

where 



.(k) 



= (i?)* =^5> CT)CT '(k)(c 

a, a' k 

+I^B r m r , (10a) 
r 

= \fq^\fq^f-(j,a'{^) + &a,cr>£a , (10b) 

= ^^'exp(-ik(R r R,„)). (11) 



The calculation only requires |$) to be a one-particle 
product wave state; |$) need not be a filled Fermi sea. 

For a given set of Coulomb parameters in JJJJ the renor- 
malization factors 

q<? = g ff (nj,mr) (12) 
only depend on the local spin-orbital densities 



(13) 



and the variational parameters mr for states |T) with 
more than one electron. An explicit expression for H12|) 
has been given in Ref. 8, but it is not needed for our 
further considerations. Note that for our symmetry- 
restricted orbital basis 



n a = (n CT )* c 



(14) 



holds. 

E vax in (|10|l depends on the variational parameters mr, 
the local densities n a , and the wave function |<3?). How- 
ever, the constraints 



E'< 



(15) 



have to be respected during the minimization as we work 
for fixed n in the sub-space of normalized one-particle 
product states, ($|$) = 1. We introduce Lagrange pa- 
rameters Egp , A(7 , and A for these constraints which leads 
to the energy functional 



. A CT , A] = S var [|$),m r ,n CT ] 

-LJ2^ [n CT -($|n CT |$)] 



-LA(n — rifj) 

(7 

-£7 BP (1 -($|$)) 



(16) 



E c has now to be minimized with respect to all quantities 
|$), Esp, mr, n a , A CT , and A independently. 

First, we use the condition that l|l(j|) is extremal with 
respect to |$) and Egp- This gives us the following effec- 
tive Schrodinger equation which has to be solved in the 
sub-space of normalized states |<E>) 



(17) 
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with 



= EE^'W + 6 °,°' X °K<A;<r> ■ (18) 



The effective one-particle Hamiltonian H eS can be diag- 
onalized, 



i^ = ££(k,r)£+A 



(19) 



by introducing proper creation and annihilation opera- 
tors 

: =EWC ; Vr-E^VA*- ( 2 °) 



Note that the amplitudes F k CT r , the energies E(k, r), and 

the operators h£ , h k still depend on the parameters 
mr, n CT , and A CT . 

Solving the eigenvalue equation H17|l is only a necessary 
but not a sufficient condition for a state |$) to minimize 
the original energy expression l |10| l. In Sect. Illll we will 
use this ambiguity to define quasi-particles excitations of 
the variational Fermi-liquid ground state. 

C. Variational Fermi-liquid ground state 

In order to obtain our variational Fermi-liquid ground 
state, it appears to be the most natural choice to proceed 
with the filled Fermi sea for the effective Hamiltonian 



|$ ) = JJ ^ka- (vacuum) 

k,r;B(k,r)<B F 



(21) 



Here, the Fermi energy Ep is determined by the condition 



j"£Q(E F -E(k,r)) 



(22) 



k.r 



The corresponding eigenvalue Esp becomes 

E S p =E- B ( k > r )9(£F-£(k,r)) . (23) 

k,r 

It is difficult to prove rigorously that the state (|21|l leads 
to the global minimum of (|l(Jfl . However, |$o) is at least 
stable with respect to single-particle excitations and it is 
difficult to conceive any other state which is consistent 
with our underlying Fermi-liquid picture. 

When we insert |$o) into (|16|) we arc led to the energy 
function 

E c [mp , n a , \ a , A] = E SP + E r m r - \ a n a 

r <t 

-LA(n-J2 n *) ■ ( 24 ) 



In the variational ground state this expression is extremal 
with respect to mp, n a , X a , and A, 



_d_ 



■E L 



{x j }={x j } 



= with Xi 6 {mp, n a , X a , A} 



(25) 

The optimum values mp, n a , X a , and A define the op- 
timum values for the energies E(\t,r), the amplitudes 

F ko . r , and the operators h£ T , /i kr . Furthermore, we 
can write the variational ground-state energy as 



E™ = E c {m r> n a) \ a ,~K\ 



(26) 



The energy (|26|l depends on the particle density n 
both implicitly, mediated by the optimum values Trip (n) , 
n a (n), \ a (n), and A(n), and explicitly, due to the term 
—LA(n)n and the Fermi energy Ep = Ep(n) in E^p 
of 124p . Therefore, the (variational) chemical potential 



can be written as 
1 dE SP 



/' 



L dn 
1 

+z 



fl = T 



A 



1 dE™ 



dn 



(27) 



(28) 



E 



dx. 



{x j }={x j } 



dxj 
dn 



The sum in l|28() vanishes due to 125(1 whereas the deriva- 
tive of E$p just gives the Fermi energy Ep, 



1 dE Si 



= E ¥ 



(29) 



L dn 

Therefore, the variational chemical potential reads 

fx = Ep — A . (30) 

The strategy for the numerical minimization is not 
important for our analysis of Landau-Gutzwiller quasi- 
particles in the rest of our work. For further reference we 
give a short summary of our most efficient procedure. 

First, we note that the conditions (dE c )/(dX a ) = 
and (dE c )/(dA) = take us back to the original con- 
straints (|15|) . 



_d 

dX, 

El* 



E SP = (?v)$ 

2 



e(Ep-E(k,v)) 



k,i 

n = E n ° 



(31a) 
(31b) 



Therefore, we are left with two different sets of variational 
parameters, the 'internal' parameters mp and the 'exter- 
nal' parameters X a . Optimizing the energy with respect 
to both of these sets is time-costly, for different reasons. 
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The problem with the internal parameters is their large 
number which is of the order 2 2N (w 500 for d orbitals). 
Compared to this there are only a few, 2N, external pa- 
rameters A CT . However, each modification of one of these 
external parameters demands for momentum-space inte- 
grations according to the sums in (|31a|) and Ij23|l . We 
found these integrations to be the most time-consuming 
part of our numerical minimization. In principle, such in- 
tegrals must also be calculated whenever we change the 
parameters mp, because they determine the amplitudes 
F kar and the energies E(k,i). In order to avoid this 
large number of integrations we write E$p in l|24|) as 

^SP = ^ x/g^y^y^ e <T, ff '( k )4*a,r^k, ( r',r 
(7,(7' k,r 

+ £(? ff + A ff )£V lw | 2 , (32) 

(7 k,r 

where the prime on the sums implies E(k, r)) < Ep. 
Eqs. (|24|) and (|32|1 show that the parameters mr enter 
the energy E c in two different ways: (i), indirectly, via 
the amplitudes F kaI or the energies E(k, r) and, (ii), 
directly, via q G in (|3*2"|) and the second term in (|2"4"}l . 
This separation suggests the following numerical itera- 
tion scheme: 

1. Start with an initial guess for the parameters mr, 
e.g., their statistical values in the uncorrelated 
limit. 

2. Minimize the energy with respect to the parame- 
ters A<j while all mr are fixed. During this mini- 
mization the constraint (|31b|) must be respected. 

3. Minimize the energy with respect to the parame- 
ters mr, while the parameters X a , the amplitudes 
F kCTr , the energies E(k,r) and the wave function 
|$o) remain fixed during step [31 

4. Go back to step |2 unless the reduction of E™ be- 
comes sufficiently small. 

The above procedure represents only a rough picture of 
our numerical minimization. For example, in practice 
one finds that some of the parameters A CT play only a mi- 
nor role and, therefore, are fixed during the whole min- 
imization. However, we are not going to discuss these 
numerical details in this work because they depend on 
the specific material under investigation. 

IV. LANDAU-GUTZWILLER 
QUASI-PARTICLES 

A. Definition 

The Gutzwiller theory provides I^q), an approximate 
description of the true many-body ground state. In or- 
der to extend the variational description to quasi-particle 



excitations, we closely follow Landau's ideas. We seek 
creation and annihilation operators e pt and v p t which 
must obey the same Fermi-Dirac distribution around the 
Fermi surface as uncorrelated electrons, i.e., we postulate 

(^° c \e+ t v T>f \^° r ) _ 

(tt°|tt°) = Q{Ef ~ S(P ' t)} ' (33) 

at zero temperature. We will see below that it is actually 
possible to define operators e pt and v pt which obey 
in the whole Brillouin zone and not only around the Fermi 
surface. This implies that our variational approach docs 
not capture the damping of quasi-particles. 

1. Quasi-particles for a rigid Fermi-sea background 

First, we adopt the viewpoint of a fixed Fermi-sea back- 
ground, i.e., we assume that a quasi-particle is added 
to the iV-particle system whose variational parameters 
have been fixed by the minimization of the energy expres- 
sion (SH i or equivalently, by the conditions 1|25|) . This 
leads to the optimum one-particle product state for the 
iV-particle system 

|$o) = ]J | vacuum) , (34) 

k,r;"E(k,r)<S F 

which, in infinite dimensions, actually is the ground state 
of the effective one-particle Hamiltonian 

t cff -^S(k,r)^ k +\ r . (35) 

k,r 

The conditions (|25[) furthermore lead to optimum param- 
eters mr and by these means define an optimum correla- 
tion operator Pq. Using the one-particle operators h pt 
and h p t we can now identify 

e+ t := fcS+^G)- 1 , (36a) 

v Pit := ^^(Pg)- 1 (36b) 

as those operators which obey the quasi-particle condi- 
tion Ij33(l . Note that the inverse operator (Pq) 1 in Ij36(l 
is well defined since we expect all parameters mr to be 
finite in Fermi-liquid systems. 

Adding/removing a quasi-particle to/from the ground 
state generates the excited states 

I^G P,t) ) q p = e p +|*° G ) =? G ^ p + |$ ) , (37a) 

I*g'V = v p , t \%) = %%, t \*o) (37b) 

with fixed Fermi-liquid background. As described in 
Sect, [n] these equations constitute an explicit example 

for Landau's ideas. The Gutzwiller correlator Pq in @ 



6 



adiabatically transforms Fermi-gas eigenstates |$) into 
(approximate) eigenstates of the Fermi liquid. 

The energy of quasi-particles or quasi-holes is defined 

as 



E qp (p,t) = ±(E™( P) t)-E?*)-n 



where 



771 var 

^0 



< r ( P ,t) = (ij)^ P , 



(38) 

(39a) 
(39b) 



The ±-sign refers to quasi-particle or quasi-hole states, 
respectively. We define the quasi-particle energy (j3*5|> in 
reference to the (variational) chemical potential of the 
system, see ((271) ■ Note that the energy in ((381) is of order 
unity whereas those in l|39|l are of O(L). 

The definition of the quasi-particle states 1(37(1 applies 
to systems of arbitrary spatial dimensions. However, in 

finite dimensions, the one-particle operators h Pt t in l|34() 
and i|37f) cannot be derived from the diagonalization of 
the effective Hamiltonian (I18f) : in this case, the opera- 
tors h Pt t must be determined by a minimization of the 
variational ground-state energy with respect to the am- 
plitudes F k(rr in PI). 



2. Quasi-particles with background relaxation 

When we add a particle to the iV-particle system we 
may expect that the variational parameters will adjust 
to the presence of the additional particle. Therefore, we 
may want to work with 



I* 



(P,t)\ 



/qp 



= P G / lp +|$o) 



P .t 



4>o 



(40) 

where the TV-particle Fermi sea |$o) is defined according 
to (|2T|> . Note that the operators hp t , h t still depend on 
the parameters mr, n a , and X a for a system with N ± 1 
particles. Then, 



£ qp (p,t) = ±(£ var (p,t) - E™) - /i 



(41) 



with 

^var 



(p,t)= Mm 

mr j A<r ,A 



(H) ( P ,t) - £A(n ± — - n„ 



is the definition of the quasi-particle and quasi-hole en- 
ergy with background relaxation. 

One may wonder whether the two definitions 1(38(1 
and (|41|t will lead to different results for the energies of 
quasi-particles and quasi-holes. Fortunately, this is not 
the case in the thermodynamic limit, i.e., 



E qp (p,t) = J B qp (p,t) + e)(l/L) 



(43) 



as we will show explicitly in appendix [S] The addi- 
tion/subtraction of one particle leads to a change in the 



optimized variational parameters to order (1/L), and, in 
principle, this could result in a change of the variational 
energy E™ to order unity. However, this quantity is ex- 
tremal with respect to the variational parameters, so that 
it changes only to order (1/L) for parameter variations 
around their optimal values. Therefore, the change of the 
quasi-particle energies due to the background relaxation 
vanishes in the thermodynamic limit. 



B. Quasi-particle dispersion 

In the following we focus on _B qp (p, t) because the eval- 
uation of 1)38(1 is more involved. The energy ((42(1 is given 
by 

E™\p,t)= Min \E^[m r ,n a ,X a ,A]\ (44) 

mr iTio- ,A CT , A L 

where 

E^[m r ,n a , A CT , A] = E^ + ^ E T m T - L ^ \ a n a 



-LK(n±j -^n a ) 



(45) 



and 

4p = ±E(p, t) + E ^ *M e f - E(k, r)) . (46) 

k,r 

The (±) sign in ((45(1 and ((46(1 correspond to a quasi- 
particle and quasi-hole state, respectively. 

Adding or removing a particle changes the parameters 
mr, n (J1 \ a , and A and the energies E(k, t) only by terms 
of the order {1/L) compared to their values in the N- 
particle ground state, 



X; X; 



with Xi G {mr, n a , X a , A} 



E(k,t) = E(k,t) 



SE(k,t) 



(47a) 
(47b) 



Thus we may expand l(45|) in terms of up to order 

unity, 



(42) ^K.s^i.A] = £ c [m r ,n a ,A„,A] 



±(E(p,t) - A) + SE C 



Xi 6 {mr ,n CT , \ a , A} 



(48) 
Sx,. 



(49) 



The sum in ((49(1 vanishes according to 1(25(1 . Using 1(29(1 
the quasi-particle dispersion 1(41(1 becomes 



Bqpfot) =E{p,t)-E F 



(50) 
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This result does not come as a surprise since the Fermi 
surfaces, as defined by the conditions E(p, t) = Ep and 
E qp (p,t) = 0, must coincide. In addition, eq. (fBTl)) shows 
that the quasi-particle dispersion is given by the eigen- 
values E(p,t) of the effective Hamiltonian l|35l) not only 
around the Fermi surface but in the whole Brillouin zone. 
Note that the variational kinetic energy, 



(Hi 



= (H 



off 



'07*0 



(51) 



is given by the expectation value of the effective Hamil- 
tonian l|35|) only in the case of degenerate orbitals where 
A<7 = 0. 

There are two important differences between the ef- 
fective Hamiltonian or, equivalently, l|18|). and the 
bare one-particle Hamiltonian Hi in (J5J. First, the bands 

are narrowed in H eS because the Coulomb interaction re- 
duces the mobility of the electrons. Second, the fields \ a 
which were originally introduced as auxiliary Lagrange 
parameters act as observable shifts of the energy bands, 
e.g., in terms of a magnetic exchange splitting. Our de- 
tailed numerical investigations on Nickel 3 ^ showed that 
both effects, i.e., band-narrowing and band-shifts, are rel- 
evant for a proper description of quasi-particles in real 
materials. 



1/D CORRECTIONS FOR THE ONE-BAND 
HUBBARD MODEL 



The energy expression (|1U|I for the wave function J2J is 
exact in the limit of infinite spatial dimensions D. There- 
fore, its evaluation for real, finite-dimensional systems 
constitutes an additional approximation. For the one- 
band model it has been shownA. 3 that 1/D corrections of 
ground-state properties are actually small in most cases 
of interest. An exception is the half-filled Hubbard model 
where, in infinite dimensions, the Gutzwiller theory pre- 
dicts the so-called Brinkman-Rice transition where all 
electrons become localized at some finite critical interac- 
tion strength Ubk- This metal- insulator is known to be 
an artifact of the limit D — > 00 because it is absent in all 
finite dimensions^. Consequently, 1 / D corrections must 
become important in this special case. 



A. First order corrections: analytical results 

In the case of only one orbital per lattice site, the gen- 
eral Hamiltonian JSJ reduces to 



# = E E e ( k )#, A, CT + u E Wf ■ (52) 



We consider a hyper-cubic lattice with only nearest- 
neighbor hopping-terms where the bare dispersion in i|52|) 



is given by 



e( k ) = ~y 7> 2^cos(fc z 



(53) 



1=1 



The Gutzwiller wave function will be evaluated in its 
original formic, i.e., the variational parameter A(-m for 
the doubly occupied state |T) = is replaced by the 
parameter g. For the one-band model both definitions 
are equivalent. 

The variational ground-state energy of the Hamilto- 
nian (1521) in infinite dimension reads 



E°°(g,n)=L[q(g,n)e + Ud(g,n)] , (54) 



where 



:= j E n k,^( k ) 



(55) 



is the mean kinetic energy of the non-interacting system. 
Here, the renormalization factors q(g, n) and the average 
double occupancy per lattice site d(g, n) are given by 

4 / Ti \ 



d(g,n) 



x[\Jl-n + d{g,n) + \Jd{g,n)) ,(56a) 

nG +n-l 
2 G+l 



(56b) 



with 



G= y/l + n(2-n)(g*-l) 



(57) 

The momentum distribution of the non-interacting sys- 
tem, 



< >CT = e(E F - e (k)) , 
and the electron density, 



(58) 
(59) 



k.a 



determine the Fermi energy E-p. 

We set up the 1/D expansion of a function A(g,n) in 
the form 



A{g,n) = A°°{g,n) + ^A^(g,n) + 



D 



(60) 



Then, the first-order correction of the ground-state en- 
ergy reads 

E^(g,n) = L[t^(g,n) + Ud (1) (g,n)] , (61) 

where the corrections to the average kinetic energy and 
the double occupancy can be written as 



k,(T 

(n-l)(G-l) 



s + e CT (k) 



(62) 



(e ) s 



n(2-n)G 
x [n(G + 1 - n) + 2(1 - n)(G - l)n Q k J , 

'(63) 
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and 

d {1 \g,n) = h(g,n)(e ) 4 
Here, we introduced the factors 



(64) 



f(9,n) 



h(g,n) = 



1 



G- 1 



1 



1 + gJ \G+1J \n(2-n) 
{G+l-n)(G + n- 1)(G- 1) 



(65a) 
(65b) 



2G(G+l) 3 n 2 (2-n) 2 
The total ground-state energy to first order in 1/D, 

E(g,n)=E°°(g,n) + ±EW(g,n) , (66) 

has to be minimized with respect to g. However, the 
optimum value g of this minimization differs from the 
respective value in infinite dimensions only by terms 
of the order 1/D, 



g = g 



l 

D 



(i) 



(67) 



Therefore, we can expand the optimum ground-state en- 
ergy in terms of 1/D as 



E(g,n) = E°°@*> + -. 

D yy 



1 



g [1 \n) 



E°°(g°°,n) 

+-U (1) or 



D' 

gW dE°°(g,n) 

~7f 



dg 



g=a" 



D 



(68b) 



which leads to 



E{g,n) » E°°(g°°,n) + ±EU(g°°,n) = E(g°° ,n) (69) 

because the derivative in (|68b|) vanishes. From (|69|l we 
see that the optimum ground-state energy is determined 
by the optimum parameter ~g°° in infinite dimensions and 
no minimization of the total energy l|66f) is required. 

In order to determine the quasi-particle dispersion as 
defined in (|38[1 . we evaluate l|39a[l and (|39b|) to order 1/D. 
The energy (|39a|) is given by l|69|) . The expression l|69|) 
also yields the energy i|39bfl when we perform the replace- 
ments 



k.tr 



-4 

77° 



(70a) 
(70b) 



for a quasi-particle state (+ sign) or quasi-hole state 
(— sign) with momentum p and spin r. A straightfor- 
ward expansion of (|39bjl in terms of 1/L leads to the 
quasi-particle energy 



E. 



qpvP.r) =E qp (p ) r) + -^(p,r) . (71) 



Here, we recover the quasi-particle dispersion in infinite 
dimensions, 

£~(p,r) - ±q(g°°,n)(s(p) - E F ) , (72) 

as already derived in Sec. IIVI The first-order correction 
reads 

E^(p,r) =±[e(p) - #f]%(p)] (73) 

(n - 1) 2 (G°° - l) 2 ^ 
n(2 - n)G° 



with 



[ £ (p)]=/(r°,n)(£o) 2 { 
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+6n(G°° + 1 - n) + 6(1 - n)(G°° - lje 2 

+ (£(p)+i? F )2(l-n)(G 00 -l)eo} (74) 

+4C/M5 oo ^)(e ) 3 . 

The quantity G°° is given by (fBTj) . evaluated at g = 7?°°, 
and E 2 , is defined by 



(75) 



Note that in deriving (|?4*|) we have used the relations 

5>(k)=0 , I]T[ £ (k)] 2 = l, (76) 

k k 

which hold for the dispersion relation (55} ■ 

B. First-order corrections: numerical results 

We are interested in the relative size of the 1/D cor- 
rections compared to the result in D = oo dimensions. 
For this purpose we introduce 



E qp p,r) 
m(p) = - 1 

^qp(P^) 



(77) 



as a measure for the deviations from the result in infinite 
dimensions. 

In the half-filled case, n — 1, the ratio m(p) is inde- 
pendent of the wave vector p, 



m(p) 



i T 



l 



2DT 



e U 



3(5°° - 1) 



r 



(78) 



The inset of Fig. ^ shows m = m(p) as a function of 
s = Ad/n 2 for spatial dimensions D = 1,2,3. Here, 
< s < 1 provides a measure for the correlation strength 
in the system. The value s = 1 corresponds to U = 
and s = is realized at the Brinkman-Rice transition, 
U = Ubr = 8|eq|- As seen from the inset of Fig. ^ 
1/D corrections are not negligible over a wide range of 



0.2 0.4 0.6 0.8 1 

S 

FIG. 1: Renormalization factor m(pp) for the quasi-particle 
dispersion at the Fermi energy as a function of s — 4d/n 2 
for band fillings n = 1, 0.99, 0.95, 0.9, 0.8, 0.5, 0.2 (from top to 
bottom at s = 0.2) in three dimensions. Inset: special case of 
half band- filling for dimensions D = 1 (dashed line), D = 2 
(dotted line), and D — 3 (full line). 

interactions, especially in one dimension. In three di- 
mensions, these corrections are much smaller but still 
about 25% close to the Brinkman-Rice transition. 

Since this transition is spurious in finite dimensions, 
1/D corrections have to be large in the half-filled Hub- 
bard model. For an application of our method to real 
Fermi-liquid systems it is more reasonable to study 
cases of non-integer band filling. Fig. ^ shows the ra- 
tio m(pp) at the Fermi surface for different band fill- 
ings n = 1, 0.99, 0.95, 0.9, 0.8, 0.5, 0.2 in three dimensions. 
The respective results for band-fillings n! = 2 — n follow 
identically due to particlc-holc symmetry. As expected, 
the corrections in Fig. ^become much smaller away from 
integer filling. 

The data in Fig. ^ show m(pp) at the Fermi energy. 
However, for n / 1, there also is a momentum depen- 
dence of m(p) which can become significant close to half 
band filling. In Fig. [21 we show the width of first order 
contributions, 

Am = Max|m(p) - m(p F )| , (79) 
p 

on a logarithmic scale for the same band fillings as in 
Fig. n Although Am strictly vanishes for n = 1 we 
see from Fig. |2 that Am can become relatively large for 
n < 1. This means that around the half-filled case we 
find 1/D corrections which strongly depend on the wave 
vector. 

As long as Am <C m(pp), a finite value of m(p) 
amounts to a rescaling of the overall band width. When 
we apply our theory to real materials^ the band width 
is basically controlled by the Racah-parameter A which 
we adjust to fit the experimental band width. There- 
fore, 1/D corrections without a significant momentum 



0.2 0.4 0.6 0.8 1 

s 

FIG. 2: Maximal width of the renormalization factor Am for 
the quasi-particle dispersion as a function of s — 4d/n 2 for 
band fillings n = 1, 0.99, 0.95, 0.9, 0.8, 0.5, 0.2 in three dimen- 
sions, shown on a logarithmic scale; notation as in Fig. ^ 

dependence will not modify the band structure in our 
variational approach. 

As shown in Figs. ^ [21 the results for D — > oo be- 
come questionable only close to integer filling and for 
very strong correlations. Therefore, we have reasons to 
believe that the quasi-particle dispersions as calculated 
in D = oo in Sect. llVl provide a good starting point for 
a sensible comparison with experimental data. The good 
agreement between experiments and our theoretical re- 
sults on nickel 3 - 4 supports such an optimistic point of 
view. 



VI. CONCLUSIONS 

In this work we used Landau's Fermi-liquid picture to 
define quasi-particle excitations in terms of Gutzwiller- 
correlated wave functions. Starting from the opti- 
mum variational ground state of a general multi-band 
Hubbard-model we identified operators which describe 
the creation and annihilation of quasi-particles in this 
state. We calculated the quasi-particle dispersion analyt- 
ically in the limit of infinite dimensions. Our variational 
states provide an illustrative example for Landau quasi- 
particles. They are also suitable for numerical investiga- 
tions, e.g., with variational Monte-Carlo techniques. 

We gave two definitions of quasi-particle operators, 
with and without a relaxation of the Fermi-sea back- 
ground. It turns out that it is more convenient to allow a 
(small) change of the variational parameters of the Fermi- 
sea background in the presence of the quasi-particles. We 
showed that both cases lead to the same result for the 
quasi-particle dispersion. This absence of a orthogonal- 
ity catastrophe is characteristic for Fermi liquids. 

Our results confirm our earlier calculations in which 
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the quasi-particle dispersion had been extracted phe- 
nomenologically from an effective one-particle Hamilto- 
nianiii£. In contrast to density-functional theory, our 
quasi-particle dispersions correspond to mathematically 
well-defined (variational) states in realistic multi-band 
Hubbard models. In general, our quasi-particle bands are 
narrower than the DFT bands because of the hopping- 
reduction factors q a in (fTTTfl . Moreover, as seen in (jT%|l . 
the Gutzwiller theory has the flexibility for the adjust- 
ment of the orbital energies through the parameters X a 
so that the DFT bands are shifted and mixed into the 
Landau-Gutzwiller quasi-particle bands. 

Our derivation of the Gutzwiller theory uses approxi- 
mations which become exact in the limit of infinite spa- 
tial dimensions, D — *■ oo. For this reason, we calculated 
first-order corrections in 1/D for the quasi-particle dis- 
persion of the one-band Hubbard model. Apart from 
the special case close to half band-filling, these correc- 
tions were found to be relatively small. Consequently, 
the quasi-particle bands as derived in D = oo for multi- 
band Hubbard models contain the essential information 
of the Gutzwiller states in three dimensions, and are thus 
suitable for a meaningful comparison with real, three- 
dimensional Fermi liquids. 
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APPENDIX A: QUASI-PARTICLE DISPERSION 
FOR A RIGID FERMI-SEA BACKGROUND 



In this appendix we evaluate the quasi-particle disper- 
sion (|38() and thereby prove that it is identical to the 
energy JSUJl. 



where 



The variational ground state 



= P G |$o) 



is given as 



(Al) 



where |<i>o) is the state (12 H evaluated for the optimum 
values mp, n a , Aer, and A. The variational ground-state 
energy l|39a|) therefore reads 



pvar 
^0 



k;o- C k;<T')$ 



(A4) 



The densities n a and the parameters rnr for the state 
(P,t)\ 



$ ' J differ from those of the iV-particle ground state 
only by terms of the order 1/L, 



1 . 

n a = n a + — dn CT , 

Li 

_ i . 

mr = mr + —omy 



where, for example, 



p,(7,t I 



(A5a) 
(A5b) 

(A5c) 



Here, the signs ± refer to a quasi-particle or quasi-hole 
state, respectively. Using (IA4|) - (|A5|) we can write the 
energy Ij39b|l as 



i?o ar (p,t) = ^ ^ Sa t g.' (k; fig + —6n a , m r + — 5m T ) 

cr.cr' k 



-L y~^ j?r(mr + ySmr) 



(A6) 



For the expectation value in (|A6() we find 

(A7) 

An expansion of (|A6|) in terms of 1 / L up to and including 
terms of order unity leads to 

S7(p,t) = E^±Y,S a Av\n,,mr)Fl^ t F^ al ^ 



<T C k;<T'/$o 
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x t: — S crrT r(k;n C r,mr) 



(A8) 



-k;CT'/$ 



cr.cr' k 



With the help of equations fT7)l and (fT%|) we find 



(A2) 



Here, we made it explicit that the numbers S a ^i(k) 
in H10|) depend on n a and mr- The state (13711 . which 
determines the expectation value (|39bfl can be written as 



I* 



(p.t). 



P G |$ 



<p.t). 



(A3) 



_d_ 



E SP = 



d 



di 



-H 



cir 



7 / $ 

^ y^( c ki(7 C k:g')¥ 
cr.cr' k 

d 

x - — S CT(T /(k;7v,m r ) 

OTly 

L(A 7 - A) 



(A9a) 

(A9b) 
(A9c) 
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where the third line (|A9c(l follows from l|24|) and (|25|l . In 

the same way we can show that 



Q 

Z2z2^ d k-,a'h -kit- Sv,a'(k;n a ,m r ) 



dmr> 

-L ■ E v < ■ (A10) 
Therefore, the energy difference in (|38|> becomes 

S™ r (p,t) - E™ = ±Y / (s^(p;n a: m r )T p ^ t F p ^. t 
+S a ^(\ a -A)\F p ^ t \ 2 ) . (All) 



The first two terms in this expression just give the eigen- 
values E(p,t) of the effective Hamiltonian (|35[1 . which 
leads to 



(A12) 



E Q ( P ,t)-£:r=±(£:(p,t)-A) 



Thus, by use of (|30[) . the quasi-particle dispersion (|39b|) 
finally becomes 



(A13) 



in agreement with our result for _B qp (p,t) as derived in 
Sec. lTVRl 
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